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Abstract 

In this paper, we develop an iterative scheme to enable the explicit calculation 
of an arbitrary post-Newtonian order for a relativistic body that reduces to the 
I Maclaurin spheroid in the appropriate limit. This scheme allows for an analysis 

' of the structure of the solution in the vicinity of bifurcation points along the 

Maclaurin sequence. The post-Newtonian expansion is solved explicitly to the 
I fourth order and its accuracy and convergence are studied by comparing it to 

■ highly accurate numerical results. 

m ■ 
o 

■ 1 Introduction 

1 ' 

^1 Upon the discovery of pulsars in 1968 and their identification as neutron stars, it became 

apparent that a relativistic description of rapidly rotating, compact stars was needed. 
Early work in this direction dealt with simplified models for the matter making up 
these objects. In particular, Chandrasekhar pP looked at stars of constant density and 
5^ i calculated the first post-Newtonian correction to the Maclaurin spheroids. Bardeen 

|2] reexamined this work using a modified approach and gained new insight regarding, 
foremost, the points of onset of secular, axisymmetric instability along a one parameter 
Maclaurin curve. 

Given the amount of work that has been done since then to study stars with more 
realistic equations of state, the return to a model of constant density can hardly be 
motivated by astrophysical considerations. Many other good arguments however, sug- 
gest that precisely this model deserves closer attention: For one, it allows, as we shall 
see, for the development of an iterative scheme to calculate explicitly any order of the 
post-Newtonian expansion, limited in practice only by computer algebra programs and 
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the machines running them. Furthermore, by considering an arbitrary order, one can 
study properties of the full relativistic solution and carry out Bardeen's task of testing 
conjectures "by going to higher orders in the relativistic expansion" . Finally, due to the 
fortuitous circumstance of being in possession of highly accurate numerical values, one 
can go even further. We are in the rare position of being able to examine the behaviour 
of the post-Newtonian expansion itself, providing, by analogy, a testbed for the most 
widely used analytic approximation within the field of General Relativity. 

In section [21 of this paper, we motivate the method used here by briefly describing 
Bardeen's approach j2] for the first order of the expansion and explaining why modifi- 
cations are necessary when going to higher orders. Section IHl presents the line element 
and the Einstein equations to be solved for iteratively. An iterative scheme allowing for 
the explicit calculation of an arbitrary order is presented in section |3] and some proper- 
ties of the solution are discussed in section El After providing by way of example the 
explicit calculation of a few expressions and introducing various physical quantities, the 
PN approximation up to the fourth order is compared with highly accurate numerical 
results in section IHl 



2 Preliminary Remarks 

For a given mass-density, Q, the Maclaurin spheroids and the relativistic model both 
depend on two parameters.^ Since the post-Newtonian approximation describes the 
relativistic model in terms of Newtonian parameters, some convention is needed to de- 
termine which relativistic parameters are implied by the specification of the Newtonian 
ones. Bardeen argued that the "most appropriate choice" compares Newtonian and rela- 
tivistic bodies of the same rest mass Mq and angular momentum J since these quantities 
(together with Q) are coordinate independent and "play the primary role in the Hartle- 
Sharp variation principle". In this paper we take a somewhat different approach 
since our purpose is less the comparison of Newtonian and relativistic configurations, 
than the development of a method for calculating an arbitrary order of the expansion. 
Therefore we use the freedom that one has in defining the PN approximation in order 
to simplify the mathematical structure of the equations. The remaining freedom re- 
garding the choice of a constant is left unspecified as long as possible. What effect the 
specification of this constant then has, will be studied in section IU7S1 of the paper. 

At this point, a brief description of the method that Bardeen used in |2j will provide 
us with the basis for understanding the motivation for the methods used in this paper. ^ 
Up to the first order of the PN approximation, one has to determine two metric functions 

^In the Newtonian case, one of these is a mere scahng parameter. 

^The reader who is interested in studying Bardeen's paper 2| rnay benefit from the following list of 
errata. The coefficient of P2{v) in eq. (21) should read —jr- 



The left hand 



side of eq. (36) should read -^{v'^)e{^s,v) and on the right hand side j2D2W2{r]). The Pi{ri) of eq. (43) 
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from Poisson-like equations as well as the unknown boundary of the star. To solve for 
the metric functions, Bardeen used a Poisson-integral in spheroidal coordinates ^ and rj, 
which represents potentials as an expansion in terms of orthogonal polynomials in rj. An 
iterative scheme for the calculation of higher orders is only feasible if the sum over these 
polynomials terminates. The conditions for the termination of the sum are that the 
source remain a polynomial in rj and that the boundary of the star remain a constant in 
^. Neither of these conditions is met with in Bardeen's approach, which is why it is only 
appropriate up to the first order. To that order it was possible to determine the metric 
functions in an elegant way, because they can be decomposed into one piece containing 
the new (post-Newtonian) source within the old boundary and another piece containing 
the old source within the new boundary. For higher orders, such a procedure can no 
longer be used and one has to devise a modified approach. 

The approach used in this paper relies on the fact that an extended version of the 
Poisson-integral is valid for Poisson-like equations even in modified coordinates. Here 
coordinates will be introduced that are tailored to the unknown boundary of the star 
and satisfy the condition that the boundary be a constant in this coordinate. Further- 
more we require that the unknown boundary of the star when written as a function of 
the old coordinate ^ be given as a polynomial in r^, a requirement that can be shown 
to be compatible with the condition that the pressure vanish at the surface of the star. 
This requirement ensures that the sources in the Poisson-like equations remain polyno- 
mials in T). Thus we have to deal only with terminating sums to any order of the PN 
approximation and the recursive method proposed here can be applied indefinitely. 

3 Basic Equations 

The line element for an axially symmetric, stationary, asymptotically flat space-time de- 
scribing a perfect fluid with purely azimuthal motion can be written in Lewis- Papapetrou 
coordinates as 

ds'' = e''^ {do" + dO + ^'e^^ {d^ - u dtf - c^e'^'dt'. 

The metric functions n, X, cu and v depend only on g and C and vanish at spatial infinity. 
The energy-momentum tensor for the pressure P and the mass density Q, which is merely 
the energy density divided by c^, is then given by 

Tap = {QC^ + P) U^UfS + QapP, 

where Q is a constant up to the surface of the star. The matter of the star rotates 
uniformly with an angular velocity fl. We introduce the spheroidal coordinates 

g' = ao'(l+^')(l~^') and C = ^0^^, ^^[-1,1], ^e[0,oo), 

is to be replaced by {Pi{r]) — 1) and the second term in the first hne of eq. (55) by +-^^(^s, 1). 
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and obtain from various combinations of the Einstein equations the following partial 
differential equations for the metric functions (with G = 1 for the gravitational constant): 



1^. 



1 + v 



r,2 



1 — V 



-iQc' + P) + 2P 



Aifi 



c4(l -{;2) 



{Qc^ + P)- C{u, 3A - u) 



—iQc' + P) + Ciu, A) + -(1 + e)il - V^''-'^^^, ^) 
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The differential operators in the above equations are defined by 

£(0, X) ■■= [(1 + + (1 - /alie + V') 

and 



(9^ /52 



d_ 



d 
dr] 



Me + 



and the dimensionless function in eq. (fig) by 



(la) 

(lb) 
(Ic) 

(Id) 



Note that the operator A2 is simply the Laplace operator in a flat three-dimensional 
space. The dimensionless pressure P := P/Qc^ is related to the metric functions by 



Vl + v^{l + P)e'' = const. = 1-7 



(2) 



with 



V := qVI{1 — uj)e^ ^/aQ and VL:=aQVL/c. 
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4 The Iterative Scheme 



4.1 The Expansion 

The system of partial differential equations is simplified by expanding the relevant 
quantities in terms of a dimensionless relativistic parameter. Here we choose the square 
root of the parameter^ used in |2] defined by 

The three variables entering into this definition completely specify the Newtonian Maclau- 
rin spheroid. Q is the mass density, ao the focus of the ellipse describing the surface of 
the star in cross-section and the value of the surface's ^ coordinate. These are the 
same quantities which will enter into the PN expansion, but the latter two lose their 
simple geometrical meaning. The parameter e remains finite in both the spherical limit, 
given by — ^ C)0 and oq oc and the disc limit, which is given by — and Q oc 
for non- vanishing mass. 

The expansion of the dimensionless metric functions and the constants reads as 
follows: 

oo oo oo 

" " " cDne" 

(4) 

As was already mentioned, Q is held constant to any order of the approximation, which 
is why it does not appear in eq. (jH) and any other quantities of interest, such as v ot P 
can be expressed in terms of these six quantities. 

If these expansions are substituted into eqs ((H), then comparing coefficients of e 
yields differential equations for the metric functions of the form A^^i = F, where 
= z/, \,Lj,fi. Because the right hand side of eqs (fTaHTcj) depends only on j > 0, 
one can solve for (pi if the lower order functions are already known. In the case of fii, 
one can calculate it from eq. ()ld|) after having determined the other three functions to 
this order, or one can compute it from an integral over 77. 

Because an analytic solution, the Maclaurin solution, for the first step is known, 
these equations would provide an iterative process for the determination of the metric 
functions to any order if the shape of the star were known. The boundary of the star 
also has to be determined iteratively however. 

•^The square root was chosen in order to enable a more convenient indexing of the expansion coeffi- 
cients. 
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We represent the surface of the star by the equation, 

(fc oo \ / oo \ 

1 + E E ^l^'(^) H ^ + E ^^(^)^' ' 
i=0 k=2 / \ k=2 / 

where we have aheady taken into account that the boundary is an elhpsoid ^ = in the 

1/2 

Newtonian order. We also require that the sum over the Legendre polynomials Cj {rj), 
a special case of the Gegenbauer polynomials discussed in section 14. 2t terminate and 
show in section U3] that this leads to a consistent solution. 



4.2 Solving the Poisson-like Equations 

In the last section an iterative scheme was proposed for the determination of the metric 
functions in which an equation of the form = F need be solved for a known function 
F = F{^,7]). The regular and asymptotically flat solution of this equation is given by 



_ _ m — 1 



1=0 



« 1 




m—1 

gnO J J KiOCi ' iv')Fi^',v')kU^\v')dv'd^' 

C -1 

In the above equation are the Gegenbauer polynomials, gj and are two linearly 
independent solutions of the (homogeneous) Gegenbauer equation defined by 



(6) 



-1 
oo 1 



with 



di' 



{l,m) ^ (0, 1) and 



(7) 



^o(0 •= arcsinh(^) 
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The term 

m 1 

kU^, f]) dr^di :=[(! + e) (1 -V^)] ' {e + V^) d^di (8) 

is a product of the volume element and the appropriate weight function for the Gegen- 
bauer polynomials and 

Kr = ^ \^ J n for m > 1 

Kl = — for / > and 



27r 



TT 



are normalizing constants. 

In eq. (jHl) the integrands jump at the surface of the star because of the jump in 
the mass density. It is therefore necessary to split them into integrals over the interior 
and exterior of the star. Clearly if the surface of the star is given, as with the leading 
order, by a constant ^ = ^s, then this division is trivial. If the boundary depends on 
Tj, then matters are complicated considerably. As of the second order in the expansion, 
the ?7-integrals no longer run over the interval rj G [—1, 1] meaning that one can no 
longer make use of the orthogonality of the Gegenbauer polynomials and one is faced 
with non-terminating sums. We hence introduce new coordinates in order to circumvent 
these difficulties. 



4.3 New Coordinates 



We introduce the coordinates 



T] = 1] 



(9) 



implying that ip = is, the boundary of the star, cf. eq. (jSj). The new coordinate ip is a 
function of both rj and e and contains the unknown coefficients 5*^^. We rewrite eqs 
in terms of the new coordinates and manipulate them such that the left hand side has 
the same form as beforehand, but with ip replacing ^. For example, the equation for 
i> = rf) now reads 



9V 



du 



(1 + ^P')—- + (1 - r/2)_ + 2^— - 2r7 



du 



These new field equations are again expanded*^ in terms of e in order to obtain a system 
of equations for 0j as was explained in section 14.11 Since ifj = ^ + 0{e'^) holds, the 

'*Note that the coefficients (j)i{tp,r]) i^ot mere transformations of if'ii^^v) since depends on e. 
Thus one must substitute ip of eq. Q into X]j=2 'Pj '7) expand the result in terms of e in order 
to express i^fc, fc < n, in terms of ^. 
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new equations for (pilip,-)]) also depend only on known functions, thereby enabling their 
recursive determination. 

The derivation of eq. (jH)) relies on the fact that in the coordinates {C,,f]), the line 
{0,1]) is identical to the line (0, —r]) and that at spatial infinity we have ^ oo. These 
properties hold for ijj as well and an analysis of the derivation shows that we are free to 
use eq. (0) as it stands, only replacing ^ by ip- 

In changing coordinates we have mapped the star onto the rectangle [0,^s] x [—1, 1], 
which means that the division of the integrals into inner and outer domains is trivial. 
The price that one pays for the simplicity in the structure of the integrals is that the 
sources of the Poisson-like equations become quite unwieldy. But the exchange of a 
conceptual for a mechanical difficulty can be termed a good deal, and all the more so 
when its result is the facilitation of the whole scheme. 



4.4 Determining the Shape of the Star 

Due to the factor in gtt of the line element, it is necessary to determine the function 
z/j_|_2 in order to calculate the i^^ order of the PN approximation. This is the only metric 
function that depends on the unknown coefficients Sji of the star's boundary.^ To 
determine these coefficients, one calculates the pressure from eq. Q along the boundary 
of the star and sets the coefficients of an expansion in terms of rj equal to zero. In 
discussing the boundary, however, it turns out to be useful to leave a portion of the 
Poisson-integral for z/j+2 unevaluated in order to arrive at the integral equation 

1 

oo „ 
1=0 

The function Bi is defined in eq. © and contains the entire dependence on Sji. The 
function bi{7]) is short for the remaining terms that result from eq. and is a known 
polynomial of order i + 2. The function fi{fj) is given by 



fiiv) ■■=9ti^s) / hfi^p) 2#(z>°)' - ^(1 - fi')iu°,y' + 2i)l - l{l + l)^z>° 



dip 



I gfW 2#(z>^)' - ^(1 - f )(z>^)" + 21)', + -^£= - l{l + l)^z>^' 



dip, 



where a dot and prime indicate partial derivatives with respect to ip and r] respectively 
and the superscripts 'i' and 'o' refer to the regions inside and outside the star. Since 



^The other metric functions depend only on Sjk, with j < i. 
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fi{r]) is a polynomial of second order, the sum in eq. pO|) terminates for polynomial 
Bi(r)). Indeed, for the form of Bi chosen in eq. (0), one arrives at a system of i + 2 
algebraic equations for i + 3 unknowns (there are i + 1 Sij to determine as well as Qi+i 
and We choose to use this system to determine all these constants but for 7i+2. 

As mentioned in section |2l this last constant can be chosen arbitrarily, which amounts 
to specifying "which" PN-approximation one wishes to have, i.e. which relativistic body 
is to be associated with a given Maclaurin spheroid. The choice of 7^+2 will be discussed 
further in section 16131 

We have shown that the form chosen for the surface of the star is consistent with 
the Einstein equations to any order of the PN expansion. This is not to say that this 
choice is unique. One can easily see that the form chosen in [2j is incompatible with 
that chosen here, since it is not a polynomial in rj. There the surface was derived having 
stipulated that the 'generating' Maclaurin spheroid should have the same rest mass 
and angular momentum as the PN star, a condition that cannot be satisfied with the 
approach chosen here. In lieu of the freedom to choose two constants, we have chosen a 
form for the boundary of the star most appropriate to our goal of devising an iterative 
scheme and can choose only one further constant. 

5 Properties of the Solution 

5.1 Reflectional Symmetry 

In Newtonian physics, it is known, that stationary, axisymmetric bodies are necessarily 
symmetric with respect to a reflection through the C = plane (see e.g. j3]). Although 
authors (e.g. [S]) have speculated that the same holds in General Relativity, it has not 
yet been proved true. In the case considered here, this symmetry arises automatically. A 
function / exhibits reflectional symmetry in (or ip-r]) coordinates precisely when it is 
an even function of t]. Because of the orthogonality of the Gegenbauer polynomials, the 
terms in the sum of eq. (jH)) for odd / are zero if F is a polynomial in 77^, a condition which 
turns out to be satisfied. The odd terms, which are provided by the unknown boundary 
coefficients Su must be zero for the boundary condition to be fulfilled. Thus we have 
shown that any axially symmetric, stationary, relativistic solution that is continuously 
connected to the Maclaurin spheroids is symmetric with respect to refiections through 
the C = plane. 

5.2 Powers of the Relativistic Parameter 

Consideration of the field equations together with the knowledge of the Newtonian be- 
haviour of the dimensionless metric functions, shows that their expansion coefficients 

^We shall see in section l5m that Sij — for odd j. 
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^21 




2 


0.17383011 


0.98522554 


0.17126187 


3 


0.11230482 


0.99375285 


0.11160323 


4 


0.08303471 


0.99657034 


0.08274493 


5 


0.06588682 


0.99783651 


0.06574427 



Table 1: Numerical values for ^^'^ the corresponding Newtonian eccentricities and 
ratios of polar to equatorial radii given by e = 1/ a/1 + and rp/r^ = ^g/ a/1 + ^s^- 



(pi begin with i = 2 and are non-zero only for even i. The same holds naturally for 7^, 
whereas Clj begins with j = 1 and appears only with odd powers powers of j. Because of 
the choice to work with the dimensionless functions introduced here, it is most appropri- 
ate to refer to the n*^ order of the PN approximation and not the half orders in between. 
What we mean by the n^^ order is that the quantities A, a), fj, and are expanded up 
to and including the order Q up to and u and 7 to 

5.3 Singularities in Parameter Space 

By comparing the highest coefficient of rj in eq. (fTTH) one can arrive at the equation 

S^^ = (11) 



with 



1 



-1 



and 



l,«+2 

n=Q 



and where /j+2 and bi are defined in eq. ()10|) . It can be shown that the denominator of 
eq. (fTT|l is proportional to the expression 

gUi^hl^iQ - es(l - esarccot(es)) =: G,{Q. (12) 

For a given (even) i, this expression vanishes for precisely one value of C,s, let us say for 
= C+2- These values, beginning with i = 2, are the points of onset of axisymmetric, 
secular instability and the bifurcation points of new axisymmetric solutions, see [3 El Ej , 
and numerical values for the first few of them can be found in Table ^ Since of 
eq. (fTTj) is not zero at the point = ^*^2) these bifurcation points are singularities 
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in the two dimensional parameter space (^s^c^o) or (^s,^)- For values of differing 
only slightly from ^*^2' configurations have properties similar to those of the 

Newtonian configurations that branch off from the Maclaurin sequence at these points. 
The Maclaurin configuration itself cannot be reached for bodies with non-zero mass 
however, and even neighbouring configurations have strict mass limitations, since e must 
be made very small in order that the PN series converge. This mass limitation can be 
inferred, for example, by referring to the tables in Appendix El Because the n^'^ PN 
order possesses a pole of order 2n — 1 at the point C, = ,^4, we expect the coefficients 
for the expansion to grow large in the vicinity of this point. ^ This is indeed the case as 
can be seen in these tables by referring to the row with = 0.17. The series containing 
these coefficients converge only for sufficiently small e as indicated above. 



6 Explicit Solution to the Fourth Order 
6.1 The Metric Functions and the Constants 

Using the iterative scheme described above, the four metric functions and the constant Q 
were explicitly solved up to the fourth post-Newtonian order. These calculations could 
in principle be carried out ad infinitum, but the lengthiness of the expressions (the fourth 
order functions would fill several hundred pages) puts a practical limit on the order that 
can be determined. Here we will merely carry out, by way of example, the calculation 
of the first few terms. 

The expansion of eq. (fTaf) with respect to the relativistic parameter e yields the 
Newtonian equations 



47r 



^2^2 = = . (13a) 



for the interior region < ^s) and 

A2Z/° = 0. (13b) 

for the region > exterior to the body. These equations are solved using eq. ^ to 
''A lengthier discussion regarding the order of the poles at ^2i+2i * > 1 can be found in [H]. 
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obtain 



+ 



(14a) 



and 



hlWCi/'iv)-hlii;)C',^\r^) 



1/2, 



(14b) 



(see Appendix 1X1 for a list of the first few gf^ and /i™). One can verify that 1^2(^3, v) = 
z/2 (^s; ''?) holds. The requirement that the pressure vanish to this order of the expansion 
then fixes the two remaining constants: 



2vi + e 



= hl{Q and 72 



(15) 



Expanding eqs (fTb|l and (fTdj) shows that X2 = f^2 = —J^2 holds and one need only expand 
eq. (fTc|l to obtain the last of the functions 02, where = z/, A, u, fi. This expansion leads 
to the equations 

A4Ci)2 = - 



-^0 Ss 



1 + 



and 



with the solutions 



A^Co^ = 0, 



t^2 



6es(i + es^)'/^^(es) + 3 (es^-a]cf(r^) 



+ 



+ ((5es^ + 4)^2'(0 + 6)]cf(ry) 



(16a) 



12 



and 



^2 -I (1 + ^s'Y^' {htio cf{rj) - hm cnri)) . 



(16b) 



Using the scheme proposed here, the calculation of the higher orders is much lengthier 
than the calculations just shown, but otherwise identical. For for example, the last 
metric function needed in describing the first post-Newtonian correction, we find 



3 (-1 - 9 + 63 77^ -0^ + 6 77^ - 54 V'^ 77^ + 3 V^) ^TT^ arccot(4) 
4 



- - (63 + 42 r;^ - 6 77^ - 9 - 54 rj' il^' + 6 77' ^ 

+ 2 772 - 36 ^2 ^2 _ ^^2 + 3 ^2 ^^2 ^ 2 ^2) ^ j^^^2 ^ ^2) es\/l+Cs') 



5. 



22 



+ 



3 (377^-1) Vl + es'arccot(es) ^ 3 277^ + 377^4' + 2 - 6 

"2 ^72 + ^^2 2 (^2 + ^2)^^^/f^ 



'02 



21 (l + V')(-l + r?)(r/ + l)^^i' 9 
4 



72 



45 (1 + 3 ^2 77^ - ^2 + 77^) arccot(^s) 

45 is' + Cs^ - + 3 77^ Cs^ + 2 77^ 
16 es'(l + es') 



(17a) 



and 



3 (-1 - 9 7/"^ + 63 7; + 6 7/2 - 54 ^ 3 V'^) \/l + Cs^ arccot(^) 

4 7/2 _|_ ^2 

3 

- -V'(27 77^ + 96 77^ -0^ + 63 77^ -0^ - 16 77^ - 78 ^2 _ 54 ^2 ^4 



+ 2 - 3 + 3 ^^) 0^ /((r?^ + ^^) (1 + ^^)^) 



5. 



22 



+ 



3 (3772 - 1) + is' arccot(^) 

2 7/2 + z/>2 



+ 



3 V; (5 7/2 + 3 ^2 ^2 _ ^2 _ 3) ^ 
2 (772 + V^2) (1 + ^2)2 



'02- 



(17b) 
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These source terms are fourth order polynomials in rj after having been multiplied with 
the factor {ip^ + rf) from eq. Because of the orthogonality of the Gegenbauer 

polynomials, z/4 is thus also a fourth order polynomial in rj. This property propagates 
itself through the successive post-Newtonian orders such that a term 0„ is always an n**^ 
order polynomial in rj. Physically this is because the perturbative-like corrections to the 
shape of the surface, which are in the form of a finite sum of Legendre polynomials (see 
eq. (0)), give rise to a finite number of multipole moments. 

The source terms in eqs (I17|l contain two of the boundary coefficients 5*^2, which 
have to be determined by solving eq. (fTn|l . i.e. requiring that the pressure vanish on the 
boundary. A term containing S12 could also be included in eqs (fTTjl . but is found to equal 
zero by applying this boundary condition. Solving eq. (jllj) and using the abbreviation 

(3 := arccot(^s) 

one finds 

^22 = ^(1 + ^s'f^''^ (288 4 /? - 45 + 408 (3 ^s' - 54 /^^ + 1575 (3^ 

- 378 13^ + 1710 (3^ - 3150 (3 - 2370 f3 - 179 
+ 1575 + 660 / (3330 (3 - 1965 + 732 (3 ^s' 

- 357 - 5075 - 3675 ^/ + 3675 /3 + 6300 p ^J^ + 27 . 

The denominator of this expression is proportional to 6*4(^3) of eq. ()12|1 and gives rise to 
the singularity at ^4, which has already been discussed. We choose to use the remaining 
two non-trivial equations extracted from eq. ()10p to determine Sq2 and ^3. Equivalently, 
we could have determined 74 instead of ^3. In either case, the remaining constant can 
be chosen freely and does not affect the validity of the solution, but instead specifies the 
Newtonian spheroid of comparison. Because the constants 5*02 and ^3 are determined 
from a linear algebraic system of equations involving S22, they also contain a first order 
pole at the point ,^4. 

The determination of higher orders proceeds identically. One first obtains the Poisson- 
like equations by expanding eqs ([TaHTcjl and extracting the coefficients of the desired 
order in e. Next one solves these using eq. (jH)) with ip replacing ^, integrating from to 
for the interior of the star and from to 00 for the exterior. The metric function /x 
can then be most easily computed by making use of the integral 

n 

/i-A = y" (/i-A),^.rfr/', (18) 
1 

where the endpoint of integration follows from lim^^i(/i — A) = (see e.g. eqs (22) and 
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(23) in [9j. The integrand can be determined from the equations 



-(/x-A),(^ = (z/ + A),fg + v^g v^Q + A,^ - yu,^ {v + A),p 



- ('^ + ^),C 7r~2 ^.C^.f' 



(19) 



and 



^(/^-A),, = ^ [('^ + A),,, - (^^ + A),cc] + \ {i^J + A,/ 



(20) 



- (/^,^?(^^ + A),e - /i,c('^ + A),c) -. ^ (tl', 

4 an 




together with the transformation equation 



/,, = aoac-«o'r/(l + a — • 



Because /x — A is a polynomial in rj to any order of the approximation, integrating eq. ()18|1 
is trivial. Finally one uses the boundary condition to determine the boundary coefficients 
and VLi+i. 

We carried out this procedure for the first four orders of the PN approximation. The 
results are entirely expressible in terms of elementary functions and in the interior of 
the star the metric functions are simply polynomials with respect to if)"^ and rf. The 
validity of the results was ensured by confirming that the disc limit of the expansion 
reduces to that of ^Hj, by showing that the expressions for the gravitational mass and 
angular momentum found by integrating over the interior of the star are identical to 
those taken from the far field and by comparing PN-values to those returned by highly 
accurate numerical calculations. 

6.2 Representative Physical Quantites 

Physical quantities that are of interest in characterizing a given configuration are its rest 
mass 




J J V i — f 

-1 



angular momentum 



1 ?B . 




II 



(22) 



-1 
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binding energy 

Eh =7MoC^ - 2nc/aoJ 

1 

- 4vrao^ J J Pe^^'^'^ [mB{v)/^sf + v'] Uv)/^s d^ ^^^^ 

-1 

and gravitational mass 

M = Mo - Eb/c^. (24) 

The expression [{i/^^b{v) / Cs)'^ + v"^] ^b{v)/^s dip drj in the integrals comes from applying 
the coordinate transformation in eq. (0) to the volume element (^^ + 1]^) d^ drj. As an 
alternative to eqs (j22|) and pljl . one can choose to calculate the angular momentum and 
gravitational mass from the far fields of uj and v respectively and then use eq. ()24|1 in 
order to find the binding energy. The disadvantage of the far field approach is that u;j+2 
must be known in order to find Jj, whereas u)j suffices otherwise. 

In addition to using the above quantities, we shall characterize configurations by the 
ratio of polar to equatorial radius 

r-p/re = = l)/V'l + (eB(r/ = 0))2, 

the polar red shift and "surface potential" Vq 

= e-''('/'=€-'?=i) _ 1 = e-^« - 1 = 7/(1 - 7), 

the central pressure 

Pc = P(V' = 0,r/ = l) 

as well as the angular velocity VL. 

In appendix |B1 tables providing information about the expansion of these quantities 
can be found. 

6.3 Convergence and Accuracy 

As was mentioned in the introduction, we are lucky to have at our disposal a highly 
accurate numerical code. The AKM code |^ uses a multi-domain spectral method 
to solve the Einstein equations for perfect fluids in an axially symmetric, stationary 
spacetime for some specified equation of state. The accuracy reached approaches ma- 
chine accuracy and has thus been used as a standard jTHI to ascertain the accuracy of 
other numerical codes such as Lorene/rotstar jHj, the SF codes [15] or that of KEH 
[T7j (see for further information). Due to the extremely high accuracy, we can 
use the numerically generated configurations as though they were analytic solutions. 
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Relative error for different choices of PN-expansion (i > 0) 



1 ^ LiliiCl IL-CLl Vdl LLC 
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/If- 
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X lU 
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n = 0.874 


















2.2 xlO- 


-5 


M = 0.004808 . . . 


3.2 X 10"^ 








3.6 xlO- 


-5 


8.6 xlO- 


-7 


2.7 xlO- 


-5 


Mo= 0.004936 . . . 


1.0 xlO~^ 


4.0 


xlO^ 


-7 


2.9 xlO" 


■4 


4.5 xlO^ 


-7 


4.1 xlO- 


-5 


Pc = 0.02151. . . 


1.2 xlQ-^ 


3.6 


xlO- 


-7 


7.7 xlO- 


-4 


1.1 xlO- 


-7 


2.2 xlO" 


-4 


J = 0.00002272 . . 


. 1.2 xlQ-^ 


1.8 


xlO- 


-6 


1.4 xlO" 


-3 






3.9 xlO" 


-4 


!i = 0.7659 . . . 

re 


4.3 xlO"^ 


3.1 


xlO- 


-7 


3.0 xlO" 


'5 


1.8 xlO- 


-7 







Table 2: Relative errors for various physical quantities according to the 4* PN order 
for different choices for the Newtonian spheroids of comparison. The header of columns 
2-6 show which equation is satisfied by the respective choice for the constant ji, i > 2. 
A dash indicates that this quantity was prescribed. 



which is what enables us to provide values for the relative errors of physical quantities 
for example. 

In what follows, we shall use units in which G = c = Q = 1 holds and use the term 
"Newtonian limit" to refer to the limit in which the Newtonian and relativistic theories 
agree. In this limit, Oq goes to zero while remains finite, thus resulting in the fact that 
M — ^ and e^" — ^ 1. A "Newtonian" or "Maclaurin" configuration on the other hand 
is the term we use to refer to the spheroidal figure that one obtains from the Newtonian 
theory, i.e. from z/2 together with 72 and Qi. 

In Table |2] a comparison was made of different choices for the Newtonian spheroid 
of comparison. This amounts to different choices for the constants 7^, i > 2. One can 
see in this example of a configuration near the Newtonian limit, that the choice made 
can lead to differences of a few orders of magnitude for relative errors. This surprising 
result can be seen, moreover, to hold over a large range of values for the parameter 
in Figure ^ On the left hand side of this figure, the PN approximation with O, = 0, 
i > 1 is depicted. The various orders react as they must in the vicinity of the Newtonian 
limit 1 — e^" = 0: each new order brings about a noticeable improvement. As one moves 
away from this limit however, the curves cross each other and it turns out that higher 
orders render a worse approximation than lower ones. The right hand side of the figure 
tells a very different story. Here the PN approximation with 7j+i = 0, z > 1 is depicted. 
Each additional term in the PN approximation brings about a marked improvement in 
accuracy and, moreover, the relative error is more than an order of magnitude lower 
than on the left hand side. 

Imagine for a moment that one had calculated the PN approximation presented here 
without being in possession of numerical values. Furthermore, let us imagine that one 
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Figure 1: The relative error of J versus 1 — e ° for configurations with Vt = 0.874. On the 
left the PN-expansion with fij = 0, i > 1, was used and on the right with •ji^i = 0. Here 
OPN refers to the Maclaurin solution, IPN to the first order of the PN approximation 
etc. 



had decided from the outset to prescribe f2j = 0, z > 1. Then one would have been 
able to produce the plot on the left hand side of Figure El without the numerical curve. 
It would have been natural to suppose that the PN series converges toward the correct 
solution and that the fourth order of the PN approximation almost provides the correct 
value for J even up to values for M of 0.12. Had one chosen 74+1 = instead, then one 
would have produced the right hand side of Figure |21 without the numerical curve and 
come to the same conclusions regarding the convergence of the PN approximation. In 
that case, however, one would have been correct. 

Although curves depicting relative errors of other physical quantities may look quite 
different from those for J shown in Figure ^ they also have many important aspects 
in common. The choice 7j+i = 0, z > 1, leads to much smaller relative errors than for 
0,1 = and one tends to find improvement with increasing order even far away from the 
Newtonian limit. These properties hold for a wide range of Q values and the relative 
errors tend to decrease with decreasing angular velocity so long as one does not come 
too close to a singularity in the parameter space. 

One well known technique for improving on the PN approximation is the use of 
the Pade approximant, which approximates a truncated series by a quotient of two 
polynomials and is discussed with reference to the PN approximation in ^Hj- For the 
disc limit of the solution considered here, it has been shown in ^UI that the Pade 
approximant provides a far better approximation of the analytic solution given in |19| 1201 
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Figure 2: A plot of J over M for configurations with f2 = 0.874. On the left the PN- 
expansion with f2j = 0, i > 1, was used and on the right with 7j_|_i = 0. Num refers to 
the numerical solution, OPN to the Maclaurin solution, IPN to the first order of the PN 
approximation etc. 



than the PN approximation itself. In the case of the Maclaurin spheroids this turns 
out to be true as well, especially for VL < 0.8. We see in Table El how well the Pade 
approximant with a polynomial of sixth order in the numerator and second order in the 
denominator converges to the correct solution. Most likely this technique would be even 
more effective when applied to a somewhat higher order of the approximation. Even up 
to the fourth order, the PN approximation turns out to be roughly comparable to older 
numerical codes even for highly relativistic configurations. An impressive illustration of 
its applicability in such highly relativistic regimes can be found in Fig. El In this figure, 
the meridional cross section of a configuration with a central pressure of 1 and a radius 
ratio of 0.7 is depicted. One can see that the surface predicted by the fourth order PN 
approximation is almost indistinguishable from the numerical values. 

For a more detailed comparison with numerical values and a more complete account 
regarding the derivation of the iterative scheme and the singularities in parameter space, 
the reader is referred to jHj 



7 Conclusion 

In this paper, an iterative procedure to enable the explicit calculation of any order of 
the PN approximation of the Maclaurin spheroids was devised. This was made possible 
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Figure 3: The meridional cross-section of a highly rclativistic configuration. For fixed 
values of and Gq, the outermost curve represents the Maclaurin spheroid, the innermost 
was calculated using the first order approximation and the curve between these using 
the fourth order approximation. The crosses are numerical values. 
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OPN 


IPN 


2PN 


3PN 


4PN 


Pade 


AKM 


a 


n 7 


n 7 


n 7 


n 7 


n 7 

U. ( 


n 7 


n 7 


n 


0.3 


0.3 


0.3 


0.3 


0.3 


0.3 


0.3 


M (xl0~2) 


7.94 


6.17 


6.244 


6.2490 


6.25013 


6.25055 


6.25070 


Mo (xlO-2) 


7.94 


7.60 


7.589 


7.5864 


7.58590 


7.58554 


7.58553 


Pc (xlO-1) 


1.47 


2.14 


2.434 


2.5627 


2.61841 


2.66099 


2.66064 


J (xlO-4) 


6.90 


6.34 


6.141 


6.0744 


6.05313 


6.04321 


6.04352 


^ (xlO-i) 


9.73 


9.79 


9.801 


9.8041 


9.80538 


9.80606 


9.80628 



Table 3: Values of various physical quantities according to different orders of the PN 
approximation. OPN stands for the Maclaurin solution, IPN for the first PN approxima- 
tion etc. The PN approximation with 7^ = 0, z > 2 was used and the Pade approximant 
was applied to the fourth order solution. 

by introducing coordinates tailored to the unknown surface of the star, by requiring 
that this surface's representation be a terminating sum and by realizing that eq. (jH)) can 
be used in the new coordinates without alteration. The PN expansion was carried out 
explicitly to the fourth order and the resulting expressions contained only elementary 
functions. 

It was proved that the n^^ PN approximation has a first order pole at ^2n+25 onset 
of the n^^ axisymmetric, harmonic mode of secular instability. The radius of convergence 
of the series becomes zero at these points, thereby dividing the ^s-^ parameter space into 
rectangles with "impermeable" walls that accumulate about (but not at) the line = 0. 
Since the PN approximation appears to converge even in the highly relativistic regime, 
it seems likely that no quasi-stationary, axisymmetric sequence of solutions leads from 
an extended, three dimensional configuration to the disc limit - all such configurations 
would have to pass through an infinite number of such impermeable walls. 

The convergence of the PN approximation was shown to depend strongly on the 
choice of the Newtonian configuration of comparison. A poor choice can render the 
approximation useless in the relativistic regime, but a good one was shown to converge 
quite well, especially when aided by the Pade approximant. These results can be taken 
as a word of warning, reminding the researcher that the PN approximation can be very 
sensitive to alterations that may have no direct physical consequences. On the other 
hand, they also demonstrate that in the best of circumstances, the PN expansion can 
yield a very good approximation to highly relativistic configurations, well beyond its 
guaranteed region of validity. 
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The functions g^{i)) and 



5o=l 

_ _ 1 

~Y~ 2 

9 35 A 15 9 3 

16 ^ 16 ^ 16 ^ 16 



9l 



h-l = arccot (■?/') 

^2 = ^ + (^-^ - arccot(V') 
,2 35^'^ 55z/; /35 ,4 15 . 3\ 



hi 



231 VJ" 119^3 231 1/' 

■ \ — H 
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/ 231 . 315 . 105 ,2 5. 



^0=1 

gl= -4V'-1 

^1 = 16^^+ 12V^ + 1 

gl= - 64^^^ - 80^^ - 24^2 - 1 



^0 = 1 



2 



4 



2 -, V'^ + 



3 3 



16 

3_16^4 12^2 1 -_V;5-4V^3_^ 



5 5 5 

64 

64 V'" 80 V' 24 V'' 1 y V' + 16V' + 8V' + V' 
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9\ 



15?/^2 3 
^ 2 



4 315 ,4 , 105 ,2 , 15 

91 = — w ^ w ^ — 

4 _ 3003 6 3465 - 945 . 35 
~ ~ T6~ ^ ~ ^ ~ Te" ^ "16 



h\= - arccot ('?/') — ^ 



2 2(^2 + 1) 

5^-^ 1\ (15^/'2 + 13)V' 



^2 = ( o o arccot(^/') + 



24 (7/'2 + l^ 



+ + -cotW (315^^ + 420^2 + 113)^ 



16 8 16y 240(7/^2 + 1) 

,4 / 429 495 ,4 135 ,2 5 \ 

hi = ^6 ^2 arccotf?/') 

V 128 128 128 128 J ^ ' 

(15015 ^6 + 27335 + 14273 ^2 + 1873) ^ 
4480 (^2 + 1) 

B Tables of Various Physical Quantities 

This appendix contains tables with the numerical values for the post-Newtonian coef- 
ficients of the quantities introduced in section 16.21 for various values of ^s- In all the 
tables, we have chosen 7, = 0, z > 2, whence we find 



7 = 72^' = 4 V 1 + (arccot(es)(l + ^s') - ^s) 



Taking into account 



2 3c 2 

87r g^sv^T+e 



we would find, for example, the following values for the second PN approximation of 
configuration with = 0.5, e = 0.7: 

ao ^ 0.32346 
e^o = 1 - 7 ^ 0.63681 

n = n/ao^ — (0.54175 e + 0.26458 - 0.32261 e^) ^ 1.285 

Co 
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M ^ (2.6180 - 1.3693 + 1.7071 e^) ^ 0.07976 
Mo ^ (2.6180 - 0.58922 + 0.99562 e^) ^ 0.08692 
Pc ~ (0.87655 - 0.29553^2 - 0.216935^) Cq^ ^ 0.07111 
(1.5346 + 1.0605 £2 + 0.243116^) ~ 0.007480 
^ fti (0.44721 - 0.55974 £2 _ 0.23039 £^) ^ 0.2283. 
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Coefficients of for the dimensionless angular velocity Q 



0.00 
0.01 
0.02 
0.03 
0.04 
0.05 
0.06 
0.07 
0.08 
0.09 
0.10 
0.11 
0.12 
0.13 
0.14 
0.15 
0.16 
0.17 
0.18 
0.19 
0.20 
0.30 
0.40 
0.50 
0.60 
0.70 
0.80 
0.90 
1.00 
1.10 
1.20 
1.30 
1.40 
1.50 
1.60 
1.70 
1.80 
1.90 
2.00 



1.0854019 

1.0716302 

1.0579561 

1.0443818 

1.0309094 

1.0175410 

1.0042785 

0.99112387 

0.97807889 

0.96514535 

0.95232494 

0.93961928 

0.92702988 

0.91455822 

0.90220565 

0.88997346 

0.87786287 

0.865875 

0.85401088 

0.84227147 

0.83065764 

0.72155838 

0.62537562 

0.54174791 

0.46979598 

0.40833633 

0.35606582 

0.3116946 

0.27402709 

0.2420016 

0.21470198 

0.19135249 

0.17130433 

0.15401892 

0.13905099 

0.12603306 

0.11466173 

0.10468625 

0.09589882 



-0.95903173 

-0.92995651 

-0.90160729 

-0.8741343 

-0.84772939 

-0.82264177 

-0.7992021 

-0.77786047 

-0.75924895 

-0.74428879 

-0.73438506 

-0.73180435 

-0.74047382 

-0.76787434 

-0.83029178 

-0.97129085 

-1.3603363 

-3.9363526 

1.9317681 

0.57651868 

0.27773007 

0.09314481 

0.1993343 

0.26458306 

0.28819495 

0.28520112 

0.26835131 

0.24576778 

0.22194776 

0.19912131 

0.17824581 

0.15961339 

0.14319157 

0.1288049 

0.11622792 

0.10523057 

0.09559856 

0.08714120 

0.07969281 



-0.21316642 

-0.19734833 

-0.18278329 

-0.16911244 

-0.15613789 

-0.14417187 

-0.13496336 

-0.13435142 

-0.16091713 

-0.28117178 

-0.8296993 

-8.5365954 

5.2676596 

4.9354647 

7.3381300 

14.854009 

49.619624 

1171.2187 

-11.780334 

13.454904 

6.1769373 

-0.2608214 

-0.43302451 

-0.32260971 

-0.17835186 

-0.06646413 

0.00483724 

0.04504743 

0.06517545 

0.07339547 

0.07493751 

0.07292802 

0.06916914 

0.06467418 

0.05999884 

0.05543697 

0.05113403 

0.04715237 

0.04350861 



-0.09051296 

-0.11039829 

-0.12768605 

-0.14231975 

-0.15411481 

-0.16467618 

-0.19414756 

-0.44265195 

-5.9366901 

12.104546 

31.560268 

809.03790 

10.231313 

-55.887234 

-134.37904 

-437.85365 

-3139.3425 

-606631.48 

2884.4373 

127.48863 

69.330710 

1.6431116 

0.15006752 

-0.27140215 

-0.34929691 

-0.28550798 

-0.19391875 

-0.11763146 

-0.06422083 

-0.02967889 

-0.00829155 

0.00454162 

0.01198574 

0.01608587 

0.01812893 

0.01891680 

0.01894809 

0.01853317 

0.0178655 



-0.04918981 

-0.03104705 

-0.01297664 

0.00475476 

0.02219393 

0.01858310 

-1.0164212 

24.323369 

413.17220 

-22.734741 

-462.39599 

-58587.792 

2190.3915 

1224.6781 

3062.2576 

15041.748 

234545.37 

3.87526090 xlO« 

-326822.46 

100.77384 

524.99073 

5.7294717 

1.3582894 

0.58262038 

0.09017830 

-0.13370899 

-0.17918978 

-0.15297576 

-0.11177801 

-0.07575039 

-0.04899455 

-0.0304128 

-0.01789323 

-0.00957769 

-0.00409702 

-0.00050821 

0.00182048 

0.00330817 

0.00423273 



Table 4: Expansion coefficients of Q for various values of with 7^ = 0, i > 2. 
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Coefficients of for tlie gravitational mass M 
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-4.6386450 


1.0140861 


1.1999932 


0.88499186 


1, 


.10 


10.182949 


-5.4737714 


0.92337032 


1.0630189 


0.86508412 


1, 


.20 


12.264778 


-6.4024064 


0.86407489 


0.94327903 


0.81216639 


i, 


.30 


14.046199 


-( .4o(030D 


O.o3193oz4 


0.543zzU ( 1 


0. /ol4zc)4z 


1, 


.40 


17.358347 


-8.5892485 


0.82324479 


0.76132333 


0.69360857 


1, 


.50 


20.420352 


-9.8700213 


0.83510223 


0.69491426 


0.64259313 


1, 


.60 


23.859349 


-11.289852 


0.86539261 


0.64129035 


0.59918663 


1 


.70 


27.700470 


-12.858921 


0.91263339 


0.59809715 


0.56295176 


1, 


.80 


31.968847 


-14.587185 


0.97583498 


0.56340574 


0.53302851 


1, 


.90 


36.689613 


-16.484440 


1.0543810 


0.53567923 


0.50848981 


2, 


.00 


41.887902 


-18.560374 


1.1479349 


0.51370787 


0.48848149 



Table 5: Expansion coefficients of — ^ for various values of with — 0, i > 2. 
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Coefficients of for tlie rest mass Mq 



0. 


.00 

















0, 


.01 


0.04189209 


0.05245596 


0.06452339 


0.07673604 


0.09064350 


0, 


.02 


0.08380931 


0.10138291 


0.12467777 


0.14552079 


0.17045889 


0, 


.03 


0.1257768 


0.14697999 


0.1810503 


0.20649775 


0.24117379 


0, 


.04 


0.16781969 


0.18946928 


0.23424172 


0.25967447 


0.30478958 


0, 


.05 


0.20996311 


0.22910739 


0.28499562 


0.30546133 


0.36899426 


0, 


.06 


0.25223219 


0.26620288 


0.33455944 


0.34955489 


0.68695096 


0, 


.07 


0.29465207 


0.30114322 


0.38575957 


0.4525081 


-5.9229587 


0, 


.08 


0.33724788 


0.3344386 


0.44661698 


2.2159335 


-123.76573 


0, 


.09 


0.38004475 


0.36679666 


0.5454014 


-3.7541975 


-0.6779804 


0, 


.10 


0.42306781 


0.39925742 


0.82632859 


-11.306658 


154.00435 


0, 


.11 


0.4663422 


0.43345453 


4.2348294 


-338.55795 


24389.392 


0, 


.12 


0.50989305 


OAizlbizA 


1 0100/1/^0 

-1.8183462 


-10.679058 


-1027.1961 


0, 


.13 


0.5537455 


0.52062728 


-1.7412446 


21.991097 


-535.19617 


0, 


1 A 

.14 


0.59792467 


0.59014308 


-3.0682419 


62.144578 


-1430.4764 


0- 


.15 


0.6424557 


0.71084828 


-7.4164107 


228.15983 


-7811.5324 


0, 


.16 


0.68736372 


0.99975724 


-28.802435 


1834.1087 


-136523.76 


0, 


.17 


0.73267386 


2.7949289 


-770.92714 


398680.63 


-2.54506740 xlO' 





.18 


0.77841126 


-1.2397211 


11.483888 


-2103.2027 


244823.57 


0, 


.19 


0.82460105 


-0.29250893 


-10.086248 


-86.857421 


270.08015 


0, 


.20 


0.87126836 


-0.07698001 


-4.9219300 


-55.908647 


-364.07926 


0, 


.30 


1.3697344 


0.01672239 


0.56758844 


-2.1242466 


-9.5006993 


0, 


.40 


1.9435987 


-0.25789658 


-1 r\ c\ -I ^ 

1.0217665 


0.07622924 


-2.9245661 


0, 


.50 


2.6179939 


-0.58921661 


0.99559791 


1.0782019 


-1.3199037 


0, 


.60 


3.4180528 


-0.92033156 


0.74161054 


1.5062120 


0.04353243 


0, 


.70 


A r\ r\ r\ r\ c\ 

4.3689082 


-1.2337985 


0.40781773 


-1 F' r\ F' A ^ 

1.5054690 


r\ F" F" F" A -1 -1 

0.95554101 


0, 


.80 


5.4956927 


-1.5299374 


0.07578544 


1.2771041 


1.3384899 


0, 


.90 


6.8235392 


-1.8155767 


-0.22248284 


0.97101036 


1.3676277 


1 


.00 


8.3775804 


-2.0987311 


-0.48106395 


0.66449818 


1.2237982 


1, 


.10 


10.182949 


-2 3866409 


-0 70436865 


38738975 


1 0182639 


1, 


.20 


12.264778 


-2.6853269 


-0.89947761 


0.14629252 


0.80533928 


1, 


.30 


14.648199 


-2.9996963 


-1.0731453 


-0.06144076 


0.60709587 


1, 


.40 


17.358347 


-3.3337847 


-1.2309272 


-0.24120271 


0.43012826 


1, 


.50 


20.420352 


-3.6909893 


-1.3771251 


-0.39850746 


0.27444199 


1, 


.60 


23.859349 


-4.0742533 


-1.5149872 


-0.53811594 


0.13764981 


1 


.70 


27.700470 


-4.4862014 


-1.6469456 


-0.66387001 


0.01679751 


1, 


.80 


31.968847 


-4.9292361 


-1.7748240 


-0.7787784 


-0.09091482 


1, 


.90 


36.689613 


-5.4056047 


-1.8999983 


-0.88516836 


-0.18789822 


2, 


.00 


41.887902 


-5.9174456 


-2.0235176 


-0.98483424 


-0.27614163 



Table 6: Expansion coefficients of —5^ for various values of ^ with 7^ = 0, i > 2. 
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Coefficients of for tlie central pressure Pc 









p2 
c 


p4 


p6 
c 


p8 





.00 


n 
u 
















0, 


.01 


n 


nnnfii s"^? 

.UUUUloO / 


-0.00172000 


0.00162459 


-0.00077003 


0.00044127 


0, 


.02 


n 


00943630 


-0.00686564 


0.00710374 


-0.00491391 


0.00493335 


0, 


.03 


n 




-0.01545281 


0.01766303 


-0.01697127 


0.02537089 


0, 


.04 


n 




-0.02755907 


0.03514794 


-0.04585596 


0.09808811 


0, 


.05 


n 




-0.04334499 


0.06242121 


-0.11038522 


0.34561893 


0, 


.06 


n 




-0.06308734 


0.10409423 


-0.25635755 


1.3551282 


0, 


.07 


n 




-0.08723187 


0.16798563 


-0.63884234 


4.1589896 


0, 


.08 


n 


03564969 

. W O tJ VJ t: Zi \J Z( 


-0.11648064 


0.26842379 


-3.0575859 


120.34604 


0, 


.09 


n 


04446798 


-0.15194212 


0.43563224 


0.41855319 


20.090197 


0, 


.10 


n 


(15412437 


-0.19540392 


0.75876395 


-2.5871083 


4.3812886 


0, 


.11 


n 


0645783Q 


-0.24986232 


2.3812791 


-100.52138 


7005.1522 


0, 


.12 


n 


0757Q584 


-0.32064223 


0.88116263 


-29.802523 


-234.29858 


0, 


.13 





08774 50Q 


-0.41805092 


2.3107894 


-34.315938 


388.81769 


0, 


.14 


n 

U . 


1 003Q634 


-0.5647382 


5.0552036 


-87.445248 


1669.9741 


0, 


.15 


n 


1 1 3791 51 


-0.82158864 


12.978500 


-348.51779 


11069.532 


0, 


.16 


n 


1 976Q41 6 


-1.4259431 


52.306181 


-3108.6680 


224222.47 


0, 


.17 





1499804'^ 


-5.0951044 


1490.6987 


-762280.11 


4.84805780 xlO^ 


0, 


.18 


n 
u 




3.0709510 


-44.367825 


4384.7277 


-570869.92 


0, 


.19 


n 

u , 




1.1138325 


19.748422 


63.519803 


-3223.0377 


0, 


.20 


n 

yj. 


. lOC/<JO'4:'4: ( 


0.63815819 


11.158362 


84.961966 


91.304718 


0, 


.30 


n 

u. 




-0.11854446 


1.3921945 


5.7408680 


18.408063 


0, 


.40 


n 


(SI nQ'^91 7 


-0.28101446 


0.41114756 


2.7524237 


4.7395486 


0, 


.50 





R7R^'^zl7 


-0.29553405 


-0.21692694 


1.7597087 


3.1658578 


0, 


.60 


1 

J. , 


1 744^^10 


-0.19762251 


-0.64456449 


0.84355303 


2.7155097 


0, 


.70 


1 


^nzLzL9(sn 


-0.01447731 


-0.86930239 


0.02192268 


1.9510756 


0, 


.80 


1 

1 , 


R(S7'^(SzLn 
.OU 1 ou^u 


0.23179265 


-0.92120251 


-0.58277995 


1.0264413 


0, 


.90 


9 




0.52592119 


-0.84400036 


-0.95208863 


0.2108996 


1 


.00 


9 




0.85851446 


-0.67570587 


-1.1260297 


-0.39050992 


1, 


.10 


3, 


.1655292 


1.2242483 


-0.44338108 


-1.1551896 


-0.78090453 


1, 


.20 


3, 


.6716491 


1.6203250 


-0.16452438 


-1.0815067 


-1.0002778 


i, 


.30 


4, 


.2159215 


z.U4o4i /U 


U.io0izU64 


c\ no A noi on 
-0.9o4y6io9 


1 nni 1 Ti T 
-i.U9ii ( i / 


1, 


.40 


4, 


.7990015 


2.4990135 


0.49416461 


-0.73571739 


-1.0879378 


1, 


.50 


5, 


.4214241 


2.9810379 


0.86390299 


-0.49678082 


-1.0156282 


1, 


.60 


6, 


.0836241 


3.4916304 


1.2572490 


-0.22655363 


-0.89172864 


1 


.70 


6, 


.7859549 


4.0310276 


1.6730773 


0.06961668 


-0.72822437 


1, 


.80 


7, 


.5287040 


4.5994981 


2.1108281 


0.38832489 


-0.53328772 


1, 


.90 


8, 


.3121058 


5.1973091 


2.5702700 


0.72740591 


-0.31249925 


2, 


.00 


9, 


.1363524 


5.8247100 


3.0513580 


1.0854878 


-0.06968863 



Table 7: Expansion coefficients of fi^a for various values of with 7^ = 0, i > 2. 
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Coefficients of for angular momentum J 



c 




c 


C 


£^ 
c 


c 


F^ 


0, 


.00 

















0, 


.01 


0.00519817 


0.01064142 


0.01709874 


0.02545363 


0.03580801 


0, 


.02 


0.01452484 


0.02955774 


0.04621122 


0.06832484 


0.09446591 


0, 


.03 


0.02637109 


0.05341626 


0.08104960 


0.11945215 


0.16152228 


0, 


.04 


0.04014030 


0.08105972 


0.11900741 


0.17572839 


0.23030652 


0, 


.05 


0.05548233 


0.11191984 


0.15848091 


0.23645053 


0.2981537 


0, 


.06 


0.07216065 


0.14575859 


0.19860914 


0.30636975 


0.52287242 


0, 


.07 


0.09000348 


0.18260101 


0.23993267 


0.43348251 


-4.1934547 


0, 


.08 


0.10888064 


0.22274743 


0.28745884 


1.8835269 


-96.001495 


0, 


.09 


0.12869081 


0.26685311 


0.36460134 


-2.9069282 


-1.4043833 


0, 


.10 


0.14935376 


0.3161068 


0.60538199 


-9.8156498 


137.90767 


0, 


.11 


0.17080538 


0.37259919 


3.8855222 


-326.89627 


23595.306 


0, 


.12 


0.19299422 


0.44011835 


-2.1803366 


-9.1109666 


-1078.5552 


0, 


.13 


0.21587897 


0.52604717 


-2.3712087 


28.000940 


-642.86950 


0, 


.14 


0.23942672 


0.64663729 


-4.3174390 


82.390748 


-1873.3462 


0, 


.15 


0.26361148 


0.84557362 


-10.556837 


317.16019 


-10852.688 


0, 


.16 


0.28841317 


1.2930930 


-41.568836 


2647.8340 


-197589.40 


0, 


.17 


0.31381674 


3.9378673 


-1128.4291 


585246.47 


-3.73910890 xlO^ 


0, 


.18 


0.33981149 


-1.9051257 


7.2345391 


-2888.5244 


340738.03 


0, 


.19 


0.36639046 


-0.4890684 


-17.726915 


-178.76932 


-44.181794 


0, 


.20 


0.39355005 


-0.13577661 


-9.2222335 


-105.05367 


-807.13712 


0, 


.30 


0.69802485 


0.52311868 


-0.57132492 


-4.7290425 


-21.880756 


0, 


.40 


1.0714436 


0.79662582 


0.05656033 


-1.1626724 


-5.0786338 


0, 


.50 


1.5346441 


1.0605245 


0.24310865 


-0.220009 


-1.5315446 


0, 


.60 


2.1149817 


1.3709717 


0.30021959 


0.05678398 


-0.10285 


0, 


.70 


2.8447368 


1.7637678 


0.33048771 


0.03348169 


0.45890186 


0, 


.80 


3.7607272 


2.2682180 


0.38836636 


-0.11923423 


0.53433897 


0, 


.90 


4.9043081 


2.9119034 


0.50500367 


-0.30541288 


0.35919637 


1 


.00 


6.3214817 


3.7232237 


0.70019629 


-0.47535233 


0.07990005 


1, 


.10 


8.0630209 


4.7327036 


0.98941174 


-0.60511924 


-0.22477905 


1, 


.20 


10.184575 


5.9736140 


1.3874168 


-0.68287074 


-0.51722445 


1 

i, 




10 7/1 




1 onooi OS 


n 7ni Q701 
-u. / Uio/ yiz 


-U. / oiUz4oy 


1, 


.40 


15.815199 


9.2980323 


2.5742218 


-0.65447365 


-1.0093842 


1, 


.50 


19.460609 


11.463577 


3.3994232 


-0.53552549 


-1.1992419 


1, 


.60 


23.758798 


14.024683 


4.4064527 


-0.33688721 


-1.3485367 


1 


.70 


28.790709 


17.030343 


5.6180320 


-0.04974338 


-1.4550308 


1, 


.80 


34.642430 


20.532741 


7.0586270 


0.33587825 


-1.5158600 


1, 


.90 


41.405215 


24.587246 


8.7544002 


0.8310431 


-1.5274045 


2, 


.00 


49.175489 


29.252409 


10.733173 


1.4478377 


-1.4852922 



Table 8: Expansion coefficients of ^^g'^s^^ for various values of with 7^ = 0, i > 2. 
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Coefficients of for tfie radius ratio r-p/r^ 







£° 


e 




^4 

e 


^6 

e 


^8 

e 


0, 


.00 


U 


















0, 


,01 


n 
U 


.uuyyyyo 


-0. 


,02602989 


0.02202631 


-0.01047418 


0.00897031 


0, 


.02 


n 
U, 


.uiyyyouu 


-0. 


,05311731 


0.05132736 


-0.03924064 


0.04708946 


0, 


.03 


n 

U, 


.UzyyoDOi 


-0. 


,08135889 


0.08911485 


-0.09392366 


0.14546497 





.04 


U 


C\OC\C\dQC\ A 

.uoyyooU4 


-0. 


,11089284 


0.13723154 


-0.18817272 


0.37882528 


0, 


.05 


U, 


.U4yyo/D2 


-0. 


,14191637 


0.19853676 


-0.3488407 


0.96345183 


0, 


.06 


n 

U, 


Ar AQnoon 

.uoyoyzzy 


-0. 


,17471238 


0.27760603 


-0.63830148 


2.8459464 


0, 


.07 


U, 


.uoyozyio 


-0, 


,20969134 


0.38211536 


-1.2643965 


7.3313198 


0, 


.08 


n 

U, 


.U ( y ( 


-0. 


,24746 


0.52597493 


-4.6830624 


172.58895 


0, 


.09 


n 
U, 


.Uoyoo ( ( u 


-0. 


,28893904 


0.73832306 


0.09914477 


26.066958 


0, 


.10 


n 

U, 


.uyyouo/ z 


-0. 


,33557622 


1.1055791 


-3.4554085 


2.2631308 


0, 


.11 


U, 


1 noQ/i n/i Q 
. iUyo4U4o 


-0, 


,38975965 


2.7622215 


-101.28858 


6954.7987 


0, 


.12 




1 1 m /I coo 
.iiyi4ozz 


-0. 


,45569145 


1.2805118 


-31.530651 


-198.73051 


0, 


.13 




.izoyiozo 


-0, 


,54145696 


2.6523173 


-35.509464 


406.54300 


0, 


.14 


n 


1 QQR/1 7fi/l 
. iOo04/o4 


-0. 


,66476338 


5.1130909 


-83.609450 


1581.1100 


0, 


.15 


n 

U, 


1 /I fiQ/1 n/i 
. i4oo4U40 


-0. 


,87312659 


11.844410 


-307.65212 


9693.8945 


0, 


.16 


n 

U, 


.io / yyuo 


-1. 


,3516343 


43.829227 


-2561.5455 


183953.19 


0, 


.17 


A 

U, 


.ib/ oyo4y 


-4. 


,2188500 


1173.6956 


-598387.12 


3.80304840 xlO^ 


0, 


.18 


(J, 


.1 M iDO 


2.1418255 


-41.159648 


3581.7333 


-454958.98 


0, 


.19 


U, 


.iooooUoo 


0.61071792 


12.872215 


12.113393 


-2462.2706 


0, 


.20 


U, 


.iyoiioi4 


0.2351988 


7.1665406 


44.350811 


-72.324822 


0, 


.30 


n 

U, 


OQ'7Q/| 7Qn 


-0. 


,37867602 


0.98503149 


1.8025725 


5.8931889 


0, 


.40 


U, 


.o / ioyUbo 


-0. 


,51715828 


0.52418984 


0.79232382 


0.41160951 


0, 


.50 


U, 


A A 701 Qi^ 

.44/ zioD 


-0. 


,5597459 


0.23038891 


0.69014814 


0.07257866 


0, 


.60 


U, 


.oi44yo ( D 


-0. 


,55020825 


0.02038324 


0.53356119 


0.27891669 


0, 


.70 


U, 


r 'TO /I OO /I 

.0 /o4dzo4 


-0. 


,51417099 


-0.10821692 


0.34347031 


0.34737635 


0, 


.80 


U, 


.Oz4oyi3Ut) 


-0. 


,46755319 


-0.17463699 


0.18619835 


0.29128276 


0, 


.90 


U, 


.DDoyt)4 ( d 


-0. 


,41919642 


-0.20141601 


0.07879986 


0.20314523 


1, 


.00 


A 

U, 


. / U/ iUb/ o 


-0. 


,37342997 


-0.20570353 


0.01287290 


0.12721316 


1, 


.10 


0, 


.73994007 


-0. 


,33202671 


-0.19838512 


-0.0247471 


0.07303315 


1, 


.20 


0, 


.76822128 


-0. 


,29543502 


-0.18575195 


-0.04470894 


0.03751262 


1, 


.30 


0, 


.79262399 


-0. 


,26347544 


-0.17123161 


-0.05416609 


0.01524089 


1 


.40 


0, 


.81373347 


-0. 


,23570866 


-0.15659797 


-0.05756331 


0.00167726 


1, 


.50 


0, 


.83205029 


-0. 


,21161961 


-0.14270895 


-0.05756365 


-0.00636027 


1, 


.60 


0, 


.8479983 


-0. 


,19070359 


-0.12992949 


-0.05573405 


-0.01094845 


1, 


.70 


0, 


.86193422 


-0. 


,17250244 


-0.11836672 


-0.05298648 


-0.01340201 


1, 


.80 


0, 


.87415728 


-0. 


,15661589 


-0.10799852 


-0.04984618 


-0.01454328 


1, 


.90 


0, 


.88491822 


-0. 


,14270141 


-0.09874277 


-0.04661101 


-0.01488493 


2, 


.00 


0, 


.89442719 


-0. 


,13046906 


-0.09049408 


-0.04344527 


-0.01474681 



Table 9: Expansion coefficients of r^jr^ for various values of ^ with 7^ = 0, i > 2. 
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